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ABSTRACT 

We use current observations on the number statistics of gravitationally lensed optical arcs towards 
galaxy clusters to derive an upper limit on the cosmological mass density of the Universe. The gravita- 
tional lensing statistics due to foreground clusters combine properties of both cluster evolution, which is 
sensitive to the matter density, and volume change, which is sensitive to the cosmological constant. The 
uncertainties associated with the predicted number of lensing events, however, currently do not allow 
one to distinguish between flat and open cosmological models with and without a cosmological constant. 
Still, after accounting for known errors, and assuming that clusters in general have dark matter core 
radii of the order ~ 35 h^ 1 kpc, we find that the cosmological mass density, fi m , is less than 0.56 at the 
95% confidence. Such a dark matter core radius is consistent with cluster potentials determined recently 
by detailed numerical inversions of strong and weak lensing imaging data. If no core radius is present, 
the upper limit on Q m increases to 0.62 (95% confidence level). The estimated upper limit on f2 m is 
consistent with various cosmological probes that suggest a low matter density for the Universe. 

Subject headings: cosmology: theory — cosmology: observations — gravitational lensing 



1. INTRODUCTION 

A large number of cosmological probes now suggest that 
\ the Universe is spatially flat with a low mass density (e.g., 

Perlmutter et al. 1998; Riess et al. 1998; Lineweaver 1998; 

Guerra & Daly 1998; Bahcall & Fan 1998). In addition to 
, the mass density, gravitational lensing statistics have al- 

■ lowed limits to be placed on the cosmological constant. 
However, current limits on the cosmological constant from 
gravitational lensing arguments are only based on lensing 

' statistics due to foreground galaxies (e.g., Kochanek 1996; 
: Falco, Kochanek, Munoz 1998; Cheng & Krauss 1999; 

■ Cooray, Quashnock, Miller 1999; Cooray 1999a; Quast & 
Helbig 1999) 1 . An alternative approach is to consider lens- 

, ing statistics due to foreground galaxy clusters (e.g., Wu 
& Hammer 1993; Bartelmann et al. 1998; Cooray 1999b). 
It is well known that galaxy cluster evolution is strongly 
sensitive to the cosmological mass density of the Universe 
(e.g., Bahcall & Fan 1998; Viana & Liddle 1998). Since 
lensing statistics are sensitive to the cosmological constant, 
it is likely that the number of lensed arcs due to galaxy 
clusters can provide strong constraints on both the mass 
density and the cosmological constant. 

Since the first suggestion that lensed optical arcs can 
be used as a cosmological probe (Wu & Hammer 1993), 
several studies have addressed specific issues related to 
the statistical calculation. These include the effect of a 
cosmological constant (Wu & Mao 1996) and background 
source evolution (Hamana & Futamase 1997). The nu- 
merical works by Bartelmann et al. (1998), using simu- 
lated clusters in three cosmological models, suggested that 
current observational statistics on lensed arcs are consis- 
tent with predictions in an open Universe (Q\ = 0) with 
~ 0.3. In Cooray (1999b; hereafter C99), we extended 

We note that other techniques, such as the luminosity distance to 
et al. 1998), also allow constraints to be placed on the cosmological cc 



the predictions to general cosmologies and also predicted 
the existence of lensed radio and sub-mm towards fore- 
ground clusters. Here, we extend the calculation in C99 
by including various uncertainties in the predicted num- 
ber of lensed optical sources to study the possibility of 
obtaining limits on cosmological parameters based on the 
observed number. 

In § 2, we describe our calculation and inputs for the 
prediction. In § 3, we compare the predicted number of 
lensed arcs to the observed number and use a reliable lower 
limit on the observed number to derive an upper limit on 
the cosmological mass density of the Universe. We follow 
the conventions that the Hubble constant, Ho, is 100ft. 
km s _1 Mpc -1 , the present matter energy density in units 
of the closure density is fijuj an d the normalized cosmolog- 
ical constant is f^A- Unless otherwise noted, quoted errors 
are l-a statistical errors. 

2. GRAVITATIONAL LENSING STATISTICS 

In this section, we briefly describe our calculation and 
especially the description of foreground lensing clusters 
(§ 2.1) and background sources (§ 2.2). We also intro- 
duce a nonsingular isothermal sphere model to describe 
galaxy cluster dark matter profile, which is primarily mo- 
tivated by recent determinations of the cluster potentials 
using high performance numerical inversions of combined 
strong and weak lensing data towards a sample of galaxy 
clusters. 

2.1. Foreground Lenses 

The differential probability that a beam towards a back- 
ground source will encounter a foreground lens with a path 

>e la supernovae at high redshifts (e.g., Perlmutter et al. 1998; Riess 
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length of dzL is: 



dr = n(z L )aicns-T—dz L , 
dz L 



(1) 



where u{zl) is the number density of foreground lenses 
at redshift zi while ai ons is the lensing cross section (e.g., 
Fukugita et al. 1992). 

Using the Press-Schechter mass function (Press & 
Schcchtcr 1974; PS), the comoving number density of 
galaxy clusters, dn(M, z), at redshift z and mass (M, M + 
dM), can be written as 



2 p da{M,z) 



ttM dM a 2 (M,z) 



exp 



dn(M, z) 
dM 



(2) 



2a 2 (M,z) 



where p is the comoving background matter density, 
a 2 (M, z) is the variance of the fluctuation spectrum aver- 
aged over a mass scale M, and S c is the linear overdensity 
of a perturbation which has collapsed and virialized. Tak- 
ing an approach similar to the one presented in Viana & 
Liddle (1998), a(M, z) is written as a function of the co- 
moving radius, R, which contains mass M at the current 
epoch: 



a(R,z) = cr 8 (z) 



R 



S/i-iMpc 



where 



j(R) = (0.3r + 0.2) 



2.92 + log 10 



- 7 (iJ) 



R 



S/i-iMpc 



(3) 



(4) 



Here, T = 0.23 ±0.05 (Peacock & Dodds 1994) is the CDM 
shape parameter; our results are insensitive to its specific 
value (e.g., Viana & Liddle 1998). In order to calculate 
growth evolution as a function of redshift in various cos- 
mologies, we write <j%{z) as: 



<7 8 (Q) g(n m (z)) 

l + z g(n m (0)) ' 



(5) 



where, following Carroll, Press & Turner (1992), the 
growth suppression factor is: 



(6) 

The normalization for as comes from the local temper- 
ature function (Pen 1998): 



/ (0.53 ± 0.05) fi-°- 46 n A = o 
ff8(Uj - 1 (0.53 ±0.05) n-°- 53 n m + n A 



i . 



(7) 



In order to model the cluster lensing potential, we use 
the nonsingular singular isothermal sphere model with the 
observed velocity dispersion and an a priori determined 
value for the core radius of the dark matter potential of the 
cluster. The evidence for a core radius in the dark matter 
profile of galaxy clusters primarily comes from the exis- 
tence of gravitationally lensed arcs in the radial direction 
from the cluster center. For simple models for the cluster 



potential involving singular isothermal models, such arcs 
are located at a distance equivalent to the core radius of 
the cluster pontential profile. Also, recent numerical in- 
versions of galaxy cluster lensing potentials using Hubble 
Space Telescope and other ground based high quality im- 
ages clearly suggest the presence of a small core radius 
(Tyson, Kochanski, DeH'Antonio 1998; Ian DcIT Antonio, 
private communication). Thus, it is necessary that we con- 
sider a lensing model which allows for the possible presence 
of a core radius. Following Hinshaw & Krauss (1987), we 
consider a isothermal sphere model with a core radius and 
write the density profile as: 



2irG(r 2 + r 2 ) ' 



(8) 



where cr vel is the dark matter velocity dispersion and r c 
is the core radius of the dark matter profile of the clus- 
ter. The conventional singular isothermal sphere (SIS) is 
recovered when r c is zero. The lensing cross section for 
the nonsingular isothermal model is given by: 



DqlDls 
Dos 



f(P) 



(9) 



where Dqli Dqs and D^g are observer to lens, observer 
to source and lens to source distances. These distances 
are calculated under the filled beam approximation. In 
Eq. 9, /(/?) is a correction factor that takes into account 
the nonsingular behavior of the density profile (see, Hin- 
shaw & Krauss 1987): 



/(/?) 



2 

1 + 5 ^T 



V^(4 + /?) 



3/2 



(10) 



where (3 is the ratio of core radius to critical radius of the 
lensing potential, with the latter measured at the redshift 
of the cluster: 



r c cH (l + z L ) ( Dps \ 
' \DlsDosJ~ 



(11) 



When the SIS model is considered, (3 = and /(/?) = 1. 
For small core radii, especially for the present case involv- 
ing galaxy clusters, one can usually ignore higher order (3 
terms associated with /(/3); we consider, however, the full 
formula in deriving cosmological parameters. Finally, the 
differential optical depth for the nonsingular isothermal 
model is: 



dr = 167r d 



^vcl(M') 



x(l + z L 



dn(M',z L ) 
dM' 



fWdM'^j (12) 



, (D OL D 



LS 



\ Dqs J dz L 



cdt 



dz L , 



The total optical depth to a given background redshift, z s , 
is given by: 

l' Zs dr 

t(z s ) = / — dz L . (13) 

Jo dz L 

In order to calculate the lensing optical depth, we take 
a two step approach to relate cluster velocity dispersion 
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to its mass. We relate velocity dispersion to cluster tem- 
perature using recently updated a — T relation (Wu et al. 
1998): 



a ve i(T) = 10 



2.57±0.03 



VkeVy 



0.56±0.09 



km s 



(14) 



and then to mass using partly theoretical M — T relation 
(e.g., Barbosa et al. 1996): 



T(M,z) = (6.8±0.5)/i3 keV 



178 

M \ * 



(15) 







(l + z). 



We have allowed for an extra uncertainty in the M — T 
relation by comparing various normalizations that have 
been suggested in the literature. Also, we note that (3 
is dependent on cluster mass through velocity dispersion 
(Eq. 11). In addition to velocity dispersion, it is likely that 
cluster core radii are also dependent on individual cluster 
masses. Even though such variations have been observa- 
tionally determined for galaxies (e.g., Lauer 1985), there is 
still no observational evidence for a dependence of galaxy 
cluster core radii with other physical properties, such as 
the X-ray luminosity or temperature. For the purpose of 
this calculation, we take a constant value for the core ra- 
dius based on the mean value of cluster core radii from 
numerical inversions (~ 35 h" 1 kpc; Dell' Antonio private 
communication). When deriving cosmological parameters, 
we vary the exact value of the core radius to investigate 
the parameter dependences on it; As we find later, our 
limits on f2 m is weakly dependent on the core radius. 




Fig. 1. Optical depth for strong lensing for a background 
source at redshift of 1, due to foreground massive clus- 
ters. Shown are the 95% confidence ranges for both flat 
(dot-dashed lines) and open (solid lines) cosmologies with 
and without a cosmological constant. Given the large un- 
certainty associated with the optical depth and the small 



difference between flat and open cosmological models, it is 
unlikely that lensed arc statistics can be used to reliably 
place limits on the cosmological constant. 

In Fig. 1, as an illustration, we show the optical depth 
to lensing due to foreground clusters with total masses 
greater than 7.5 x 10 h~ x Mq for a background source 
at a redshift of 1, as a function of Cl m for open and flat 
cosmologies, and considering a lensing potential in which 
f(/3) — 1 (SIS model). Shown are the 95% upper and lower 
confidences in each case by considering all possible errors 
we have so far considered. The uncertainty in the optical 
depth is primarily dominated by the error associated with 
the normalization of the PS mass function; since the num- 
ber density of massive clusters is strongly sensitive to the 
exponential term in Eq. 2, small changes in erg can pro- 
duce order of magnitude changes in the number density. 
The difference between flat and open cosmological models 
is primarily due to the increase in lensing probability with 
the addition of However, this difference is small, and 
when errors in observations are also considered, it is im- 
possible to study the possible existence of a cosmological 
constant using lensed arc statistics. Therefore, taking a 
conservative approach, we combine the upper curve valid 
for flat cosmologies with the lower curve defined by open 
models to combine the 95% confidence range in the pre- 
dicted number of lensed sources. 



2.2. Background Sources 

In order to obtain reliable predictions on the number 
of lensed arcs, it is important that both the background 
source evolution and effects such as "magnification bias" 
(Kochanek 1991) be included in the calculation. Our de- 
scription of background sources comes from the Hubble 
Deep Field (HDF; Williams et al. 1996). We use the HDF 
redshift and magnitude distribution and the luminosity 
function from Sawicki et al. (1997). Such an approach al- 
lows us to reliably account for the true redshift distribution 
of background sources, instead of an empirical distribution 
or a constant redshift, while also accounting for intrinsic 
evolutionary effects which has shown to be important for 
lensing predictions (e.g., Hamana & Futamase 1997) 

Using the probability, t(z, Q m , M m ; n ), for a source at 
redshift z to be strongly lensed and the number of unlcnscd 
background sources between rest-frame luminosity L and 
L + dL and between redshifts z and z + dz, $(L_, z)dLdz, 
we can write the number of lensed galaxies, dN, in that 
luminosity and redshift interval as (see, also Maoz et al. 
1992): 



dN(L, z) 



dz 

«!.( V dL 



r(z,r2 m ,M min ) 



(16) 



A 



A 



f(A,L,z)q(A) dA. 



Here, the integral is over all allowed values of A, the am- 
plification of the brightest lensed image, q(A) is the prob- 
ability distribution of amplifications, and f(A,L,z) is the 
probability of observing the brightest image given A, L, 
and z. Our assumption that the lenses are nonsingular 
isothermal spheres implies that the minimum amplifica- 
tion, A m i n , is a function of (i. In general, the probability 
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distribution of amplifications can be written as: 
q(A) dA = 2A 2 min A- 3 dA. 



(17) 




Fig. 2. Minimum amplification versus (3, the ratio of 
core radius to critical radius at the redshift of the lensing 
cluster. 



where the sum is over each of the background galaxies. 
The index i represents each galaxy; hence, z i? Lj, and rrij 
are, respectively, the redshift, rest-frame luminosity, and 
apparent magnitude of the ith galaxy. 

Since Li for individual galaxy is unknown, due to uncer- 
tain K-corrections, following Cooray, Quashnock & Miller 
(1999), we estimate the total average bias by weighting the 
integral in Eq. 19 by a normalized distribution of luminosi- 
ties Li drawn from the Schechter function appropriate for 
the redshift Zi of galaxy i. We calculated the magnifi- 
cation bias for individual redshift intervals for which the 
Schechter function parameters are available in Table 1 of 
Sawicki et al. (1997). In principle, the uncertainties in 
the Schechter function parameters at a given redshift can 
affect the calculation of the bias, but in practice only the 
uncertainty in the power-law slope a has a significant ef- 
fect. The effect of varying a on the lensing statistics was 
discussed in Cooray, Quashnock & Miller (1999) for lens- 
ing statistics involving foreground galaxies in the Hubble 
Deep Field, which is also valid for the present case involv- 
ing galaxy clusters; the general effect due to uncertainties 
tabulated in Sawicki et al. (1997) is that the constraints 
on cosmological parameters vary by less than 5% percent, 
when a is in general varied by the quoted 1-tr uncertain- 
ties in Sawicki et al. (1997). Here, we take a conservative 
approach and allow the largest possible bias, so that the 
expected number is overestimated by an amount as sug- 
gested above. The only effect of this approach is to slightly 
increase our upper limit on fl m . 



In Fig. 2, we show A m i n as a function of (3, which is cal- 
culated following Cheng & Krauss (1999). In practice we 
use a fitting function that returns A m i n for a given value of 
(3, with an accuracy of better than 0.1% at all interested 
values of (3 in the present calculation. For simplicity, we as- 
sume that f(A, L, z) is a step function, 6[mn m , A], so that 
a lensed image with apparent magnitude brighter than 
mi; m is detected. For a given value of the core radius and 
the velocity dispersion, (3 is determined from Eq. 11. For 
massive clusters discussed here with velocity dispersions of 
the order ;> 1000 km s" 1 and at redshifts ~ 0.2, /3 £ 0.07 
and A m i n <J 5. Compared to the SIS model, the addi- 
tion of a small core radius only produces slight changes in 
the lensing probability. As described in Kochanek (1995), 
the effect of a core radius is to increase the magnification 
bias while increasing the effective lensing cross section; the 
overall effect is that the presence of a core radius is not 
significantly different from that of a SIS model. 

We assume that the brightness distribution of back- 
ground galaxies at any given redshift is described by a 
Schechter function (Schechter 1976), in which the comov- 
ing density of galaxies at redshift z and with luminosity 
between L and L + dL is 



<P{L,z)dL = <j>*(z) 



L 



L*(z) 



e -L/L' {z) dL 



(18) 
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where, as before, both L and L* are measured in the rest 
frame of the galaxy. Following CQM, we can write the 
expected number N of lensed sources as 

l-a(zi) e Li/ L' { Zi ) e -Li/AL* ( Zi ) 



10" 



Fig. 3. Expected number of lensed arcs on the whole 
sky with amplifications greater than 10 and V-band mag- 
nitudes brighter than 22 towards foreground massive clus- 
ters. The shaded range shows the 95% confidence upper 
and lower limits on the expected number of lensed arcs, 
while the horizontal lines show the range of current ob- 
served numbers. We use the lower limit on the current 
observed number to impose an upper limit on Q m . 

In Figure 3, we show the expected number of lensed 
arcs towards foreground massive clusters with total mass 
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greater than M m j n = 7.5 x 10 h Mq, and assuming 
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: dA 



(19) 
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a zero core radius for the lensing potential. We define 
an arc as a lensed source which is amplified by a factor 
equal to or greater than 10. To make a direct comparison 
to both observations and prior predictions, we impose a 
limiting V-band magnitude of 22. Our numbers can be 
directly compared to previous estimates, especially those 
of Bartelmann et al. 1998). This study predicted ~ 2400 
arcs in an open Universe with Q m <~ 0.3 and <~ 36 arcs in 
an Einstein-de Sitter Universe. The number expected in a 
flat Universe with Cl m ~ 0.3 and J7a ~ 0.7 was <~ 280. Our 
estimates for for an Einstein-de Sitter Universe range from 
~ 0.1 to 60 while for Q m ~ 0.3 Universe (independent of 
Oa) is ~ 50 to 7000 (with the higher end allowed by the 
cosmological constant). The primary reason for a lower 
number of arcs with J7a in the study by Bartelmann et 
al. (1998) was their assumption that clusters are different 
in Universes with a cosmological constant, such that their 
concentration is lower. Based on numerical simulations 
performed by the Virgo Consortium, however, Thomas et 
al. (1998) studied a series of clusters in four different cos- 
mologies, including an open model with TO = 0.3 and 
a flat model with a cosmological constant of f^A = 0.7. 
The authors concluded that clusters do not exhibit differ- 
ences between open and flat cosmologies with and without 
a cosmological constant and that cluster structures cannot 
be used to discriminate between the two possibilities. If 
Thomas et al. (1998) are correct, then the inclusion of 
a cosmological constant is not expected to change cluster 
mass profiles to an extent that would affect the gravita- 
tional lensing rate. In any case, such systematic effects are 
unlikely to be nearly as large as the current uncertainty 
in £78 which dominates the present calculation on the lens- 
ing rate. Ignoring this case, our predictions are generally 
consistent with Bartelmann et al. (1998). As we have 
demonstrated in Fig. 1, lensed arc statistics are unlikely 
to provide useful limits on the cosmological constant. The 
same is true for alternatives to the cosmological constant, 
such as scalar field and quintessence models that have re- 
cently been introduced (e.g., Stcinhardt et al. 1998). 

3. CONSTRAINTS ON n M 

In order to derive a limit on Cl m based on the number 
of lensed arcs, we require knowledge on the observed num- 
ber of such lensing events. Current surveys of clusters are 
based on their X-ray luminosities rather than masses. For 
an example, the luminosity cutoff of the followup Einstein 
Medium Sensitivity Survey (EMSS) cluster arc survey by 
Le Fevre et al. (1994) is 8 x 10 44 /i~ 2 ergss _1 , measured 
in the EMSS band of 0.3 to 3.5 keV. Converting this to a 
total mass by following through a recently derived L — T 
relation (Arnaud & Evrard 1998) in addition to the above 
M — T relation suggest a reliable lower limit on the mass of 
7.5 x 10 14 h^ 1 Mq. We have taken the lowest limit on mass 
by considering all cosmologies - since the M — T relation 
and L estimates are different under varying cosmologies. 
The current observed arc statistics (e.g., Le Fevre et al. 
1994; Luppino et al. 1998), when converted to a whole sky 
number, suggest that the number of arcs towards above 
defined massive clusters and with amplifications greater 
than 10 down to a V-band limiting magnitude of 22 is be- 
tween 1500 and 2500 (e.g., C99; Bartelmann et al. 1998). 
Ignoring the upper value, which is likely to be unreliable, 



we use the lower estimate to derive an upper limit on Q m . 
Since the lower estimate is based on the observed number, 
this allows us to put a reliable upper limit on fl m . We also 
vary this lower limit to study its effects on our constraints. 

In order to derive a constraint on Cl m , we adopt a 
Bayesian approach, and take a uniform prior for Q m be- 
tween and +1. This is primarily due to the fact that we 
do not yet have a precise determination of Q m , and, based 
on various theoretical arguments, we do not wish to con- 
sider cosmologies in which either this quantity lies outside 
the interval [0,1]. Since the prior for O m is uniform, the 
posterior probability density is simply proportional to the 
likelihood. 

The likelihood C — a function of Cl m — is the prob- 
ability of the data, given f2 m . The likelihood for n ob- 
served arcs (at redshifts Zj) when N is expected is given 
by (Cooray, Quashnock & Miller 1999): 



(C(n)) = '[[T(z j )xe-*x(l + v. 



3=0 



N 2 - n(n- 



i) 



(20) 

We have taken into account the uncertainty in predicted 
lensing rate by introducing a T , which is the fractional 1-cr 
error on t and then mariginalizing the likelihood over the 
variance of it. In order to constrain Ct m , we also need the 
redshifts Zj of the observed arcs. Using published data on 
individual lensing clusters, we obtained a median redshift 
for the lensed arcs of z ~ 1.6. As we find below, changing 
this mean redshift to a reasonably different value does not 
change our constraints on the fl m greatly 

When the observed lower limit is compared with predic- 
tions, we find that fl m <J 0.62 at the 95% confidence when 
there is no core radius (SIS). When we include a core ra- 
dius of 35 h -1 kpc, the upper limit on Q, m decreases to 0.56 
at the 95% confidence level; The change in the upper limit 
on Vt m is only a minor effect. If the core radius were to be 
as large as 100 h -1 kpc, then the derived upper limit on 
fl m can be as low as 0.29 at the 95% confidence level. How- 
ever, such a large core radius for the cluster dark matter 
profile is ruled out leaving the possibility only for a much 
smaller core radius of the order 30 to 40 h _1 kpc. When 
the effective median redshift of lensed arcs are changed to 
a lower number as <~ 1, the upper limit increases to 0.58 
from 0.56, while when the redshift is increases to a value 
of 3.0 from 1.6, the upper limit on Cl m at the 95% confi- 
dence decreases to 0.52. When we increase the lower limit 
from 1500 to 2000, our upper limit on fi TO with a model 
involving core radius of 35 h _1 kpc, decreases to 0.52 from 
0.56. This is primarily due to the fact that the expected 
number of lensing events varies by orders of magnitude 
when O m is changed from 1 to 0, with the variation in the 
expected number larger at the lower end of Q m values. For 
such a small core radius, the limit on il m is consistent with 
current estimates based on other cosmological probes such 
as type la supernovae and galaxy cluster abundances. We 
note that our limit on Q m does not mean that the Universe 
is open without a cosmological constant, but rather lensed 
arc statistics are not sensitive enough to the cosmological 
constant to see its effects above the current uncertainties. 
In general, the upper limit on Q m with a cosmological con- 
stant is slightly higher when compared to an open model. 
However, this difference is rather small (~ few percent, see 
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Fig. 1), and cannot be distinguished using current obser- 
vations on cluster number counts and lensed arcs. 

3.1. Uncertainties & Systematic Effects 

Using a lower limit on the observed number of lensed 
arcs, we have derived an upper limit on f2 m . A major un- 
certainty is likely to come when estimating a lower limit on 
the observed number of arcs since it is only based on op- 
tical followup observations of EMSS clusters (e.g., Henry 
et al. 1992). As a reliable approach, we have taken the 
lower limit allowed by the observed number of lensed arcs 
towards this sample. In reality, the true number is likely 
to be higher but the lower limit allows us to safely con- 
sider upper limits on cosmological parameters, especially 
the cosmological mass density. The present observational 
number on the number of lensing events is unlikely to be 
improved unless large samples of clusters are followed up 
at optical wavelengths. Several attempts are currently un- 
derway (e.g., Luppino et al. 1998), however, all such sur- 
veys our still based on the EMSS sample. It is likely that 
the optical followup observations of additional cluster cat- 
alogs, such as the ROSAT Bright Cluster Survey (BCS; 
Ebeling et al. 1998), can greatly improve our knowledge 
on the lensing statistics due to galaxy clusters allowing 
better constraints on the cosmological parameters. 

In addition to current low number statistics, other un- 
certainties are likely to come from the conversion of obser- 
vations, such as cluster X-ray luminosity, to mass. How- 
ever, at each step, we have considered various estimates 
such that the predicted number of lensed arcs is overesti- 
mated; This approach allows us to consider a reliable limit 
on fi m , whose upper limit may have been systematically 
increased by our procedure. We have also investigated the 



effect of a core radius on arc statistics. As found, for lu- 
minous optical arcs with amplifications greater than 10, 
the effect of a core radius on our prediction on the number 
of lensing events is minimal. The upper limit only varies 
from 0.62 to 0.56 at the 95% confidence when a reasonable 
core radius of size 35 h" 1 kpc is introduced. Increasing the 
core radius as high as 100 h _1 kpc reduces the upper limit 
by a factor of ~ 2, however, such a large core radius is 
ruled out by current observations of gravitational lensing 
of clusters (e.g., the nonexistence of radial arcs at large 
distances from the cluster center). 



4. SUMMARY & CONCLUSIONS 

Using a lower limit on the observed number of lensed 
arcs due to clusters, we have calculated an upper limit 
on fi m . Due to large uncertainties in the predicted num- 
ber of lensed sources, primarily dominated by the error in 
erg, we are unable to place limits on the cosmological con- 
stant. However, after considering possible known errors, 
and carefully taking account various estimates such that 
the upper limit on fl m is not reduced, we conclude that 
Q m <; 0.62 at the 95% confidence. 

I am grateful to Richard Mushotzky for pointing out 
the possibility to derive an upper limit on £l m , Ian 
Dell' Antonio for communicating details of his numerical 
inversions and constraints on cluster density profiles. I also 
acknowledge useful comments from an anonymous referee 
which led to several improvements in the paper and ac- 
knowledge partial support from a McCormick Fellowship 
at the University of Chicago. 
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